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We propose a multiconfigurational hybrid density-functional theory which rigorously combines a 
multiconfiguration self-consistent-field calculation with a density-functional approximation based on 
a linear decomposition of the electron-electron interaction. This gives a straightforward extension 
of the usual hybrid approximations by essentially adding a fraction A of exact static correlation 
in addition to the fraction A of exact exchange. Test calculations on the cycloaddition reactions 
of ozone with ethylene or acetylene and the dissociation of diatomic molecules with the Perdew- 
Burke-Ernzerhof (PBE) and Becke-Lee- Yang-Parr (BLYP) density functionals show that a good 
value of A is 0.25, as in the usual hybrid approximations. The results suggest that the proposed 
multiconfigurational hybrid approximations can improve over usual density-functional calculations 
for situations with strong static correlation effects. 



C/3 

Or 



I. INTRODUCTION 

Density-functional theory (DFT) [1] within the Kohn- 
Sham (KS) scheme [2] is the most widely used method 
for electronic-structure calculations in atomic, molecu- 
lar and solid-state systems. With the usual approx- 
imate density functionals, such as generalized-gradient 
approximations (GGA) and hybrid approximations mix- 
ing in a fraction of Hartree-Fock (HF) exchange, DFT 
KS generally gives good results for situations in which 
the so-called dynamic electron correlation dominates the 
total correlation energy, but it can yield severe errors 
for systems with strong static correlation, i.e. with par- 
tially filled near-degenerate orbitals (see, e.g., Ref. 3). It 
has been argued that the GGA exchange density func- 
tionals actually mimic some static correlation through 
their self-interaction error, though in an imperfect man- 
ner (see, e.g., Refs. 4, 5). There is often indeed a partial 
cancellation of errors between the self-interaction error 
which tends to give too low an energy and the neglect 
of static correlation which gives too high an energy (see, 
e.g., Ref. 6). Hybrid approximations have a smaller self- 
interaction error and are thus often worse than pure den- 
sity functionals for describing systems with static corre- 
lation (see, e.g., Ref. 7). 

Several approaches have been proposed to include ex- 
plicit static correlation in density-functional theory. Ar- 
tificially breaking (space and spin) symmetry by unre- 
stricted KS calculations is the simplest approach to simu- 
late static correlation (see, e.g., Ref. 8), and it often leads 
to reasonable potential energy surfaces but wrong spin 
densities. Another possible approach consists in replac- 
ing the single KS determinant by an ensemble of determi- 
nants or, equivalently, using fractional occupation num- 
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bers for the orbitals [9-22], but a successful and general 
method based on this idea is still lacking. Configuration- 
interaction schemes have also been proposed in which 
modified Hamiltonian matrix elements include informa- 
tion from DFT [23-26]. A lot of approaches consist in 
adding to the energy of a partially correlated wave func- 
tion calculation, including near-degenerate configuration 
state functions coupled by the full Coulombic electron- 
electron interaction, an energy density functional describ- 
ing the missing correlation effects [27-62]. In these last 
approaches, one must use a density functional which de- 
pends on the size of the multiconfigurational expansion, 
in order to avoid double counting of correlation between 
the wave function part of the calculation and the den- 
sity functional. Finally, to avoiding any double count- 
ing of correlation from the beginning, it has been pro- 
posed to decompose the Coulombic electron-electron in- 
teraction into long-range and short-range components, 
the long-range part being treated by a method capable 
of describing static correlation and the short-range part 
being described by a density functional approximation. 
The methods that have been used for the long-range part 
are: configuration interaction [63-69], multiconfigura- 
tional self-consistent field (MCSCF) [70-72], multirefer- 
ence perturbation theory [73] , constrained-pairing mean- 
field theory [74, 75], and density-matrix functional the- 
ory [76-78]. 

In this work, we explore the possibility to combine MC- 
SCF and DFT based on a simple linear decomposition of 
the Coulombic electron-electron interaction, in the spirit 
of the usual hybrid approximations, and similarly to what 
was recently done for constructing theoretically justified 
double-hybrid approximations [79]. This approach gives 
a straightforward multiconfigurational extension of the 
standard hybrid approximations, and aims at improv- 
ing their description of static correlation. After deriv- 
ing this multiconfigurational hybrid density-functional 
theory, we test this approach on situations with strong 
static correlation effects, namely the cycloaddition reac- 
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tions of ozone with ethylene or acetylene and the dissoci- 
ation of diatomic molecules, and we compare with other 
methods, in particular the range-separated multiconfigu- 
rational hybrid method of Ref. 70-72. 

II. THEORY 

Using the formalism of the multidetcrminant extension 
of the KS scheme (see, e.g., Rcfs. 69, 71, 79, 80), for any 
coupling constant A, the exact energy can be expressed as 
the following minimization over multidetcrminant wave 
functions 'J: 

E = min{(#|f + V ext + XW ee \^)+E^ :c [n 9 }}, (1) 

where T is the kinetic energy operator, V CK t is a 
scalar external potential operator (e.g., nuclei-electron), 
W ee is the electron-electron interaction operator, and 
Ey xc [n$] is the complement A-dcpcndent Hartrec- 
cxchange-correlation density functional evaluated at the 
density coming from The complement density func- 
tional, Eft xc [n} = Eu xc [n] - E^ xc [n], is the difference 
between the usual KS density functional E-R XC [n] and the 
A-dependent density functional E^ xc [n] corresponding to 
the interaction XW ee - The Hartree-cxchange contribu- 
tion is of first order in the electron-electron interaction 
and is thus linear in A, 

E^ x [n] = (l-X)E Hx [n], (2) 

where .Eh^M is the usual KS Hartree-exchange density 
functional. The correlation contribution is obtained by 
uniform coordinate scaling of the density [81-84], 

E x c [n] = E c [n]-E x c [n] 

= E c [n}-X 2 E c [n 1/x ], (3) 

where E c [n] is the usual KS correlation functional, E* [n] 
is the correlation functional corresponding to the inter- 
action XW ee , and ni/\(r) = (1/A) 3 n(r/A) is the scaled 
density. 

The theory is so far exact but in practice approxima- 
tions must be used for the multidetcrminant wave func- 
tion and the density functionals. In Ref. 79, by restricting 
the search in Eq. (1) to single-determinant wave func- 
tions <E>, we defined the density- scaled one-parameter hy- 
brid (DS1H) approximation 

£ D S1H = min (($|f + K ext + XWeeW) 

+ (1 - X)E Hx [n*) + £ c [n$] - A 2 £ c [n $ , 1/A ]}, (4) 

and, by additionally neglecting density scaling in the cor- 
relation functional, E c [ni/\] sa _E c [n], we obtained the 
one-parameter hybrid (1H) approximation, 

E m = min{(*|T + V cxt + AT^ ee |*> 

+(1 - X)E Hx [n^] + (1 - A 2 )£ C M}, (5) 



which is similar to the usual one-parameter hybrid ap- 
proximations [85, 86], except that the correlation func- 
tional is weighted by a factor of (1 — A 2 ). Starting 
from these references and applying a second-order M0ller- 
Plesset (MP2) perturbation theory [80, 87, 88], we also 
defined the density-scaled one-parameter double-hybrid 
(DS1DH) and one-parameter double-hybrid (1DH) ap- 
proximations [79], which are one-parameter versions of 
the original double- hybrid approximations [89]. These 
latter also combine HF exchange and MP2 correlation 
with a scmilocal exchange-correlation density functional 
but with two empirical parameters. 

Here, we follow a different route and use an MCSCF 
wave function in Eq. (1), expanded as a linear combina- 
tion of configuration state functions 

|*) = 5>|*/>, (6) 
i 

where the coefficients c/ and the orbitals are to be si- 
multaneously optimized. With this form of wave func- 
tion, we obtain a multiconfigurational density-scaled one- 
parameter hybrid (MCDS1H) approximation, 

^MCDSIH = min ((^|f + Fe Xt + XWee\*) 

+ (1 - X)E Hx [ny] + E c [n*} - A 2 £ C K, 1/A ]}, (7) 

and, if density scaling, which is not considered in usual 
hybrid approximations, is neglected, we obtain a multi- 
configurational one-parameter hybrid (MC1H) approxi- 
mation, 

E Mcm = min{(*|f + V ext + XW e e\^} 

+(1 - A)£ H *K] + (1 - A 2 )£ C K]}. (8) 

Equation (7) is the equivalent of the range-separated mul- 
ticonfigurational hybrid method of Rcfs. 70-72, that we 
will refer to as MC-srDFT, but for a linear separation of 
the electron-electron interaction. Notice that, if we were 
to use no approximations for the wave function 'J and 
the exchange-correlation density functional, then Eq. (7) 
would give the exact energy, independently of A. In prac- 
tice, of course, we must use approximations, and the en- 
ergy does depend on A, which can then be considered as 
an empirical parameter to be optimized. 

The present scheme has two advantages over the range- 
separated scheme: (a) Only one list of two-electron 
Coulomb integrals is needed and it is just multiplied by 
A in the MCSCF part and by (1 — A) in the complement 
Hartree energy, whereas two lists of two-electron inte- 
grals are needed in the range-separated scheme for the 
long-range MCSCF part and for the short-range comple- 
ment Hartree energy; (b) No new exchange and correla- 
tion density functionals need in principle to be developed 
since all the existing approximations developed for the 
KS scheme can be reused with a simple scaling, whereas 
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new short-range density-functional approximations must 
be developed in the range-separated scheme. 

Equations (7) and (8) can be seen as straightforward 
multiconhgurational extensions of the usual hybrid ap- 
proximations. Indeed, the expectation value of XW ee 
over the MCSCF wave function ^ not only introduces 
a fraction A of exact exchange but also a fraction of ex- 
act static correlation. Defined in the ideal limit of refer- 
ence energy levels that are degenerate, the static correla- 
tion energy is linear with respect to the electron-electron 
interaction, and thus we can consider that the expec- 
tation value of XW ee introduces a linear fraction A of 
static correlation. By contrast, the dynamic correlation 
energy starts at quadratic order in the electron-electron 
interaction, so that, for sufficiently small A, it is justi- 
fied to neglect it in the wave function expectation value. 
Moreover, for sufficiently small A, the weight (1 — A 2 ) is 
close to 1 and thus Eq. (8) includes a nearly complete 
approximate correlation energy functional, that is often 
thought of as correctly describing dynamic correlation. 
Of course, if the multiconfigurational hybrid approxima- 
tions of Eq. (7) or (8) are to be accurate, the fraction 
(1 — A) of static correlation energy not treated by MCSCF 
must be accounted for by the density functional, possibly 
through a compensation with the self-interaction error of 
the scaled exchange functional (1 — X)E x [n]. 

III. COMPUTATIONAL DETAILS 

The calculations have been performed with a devel- 
opment version of the D ALTON 2011 program [92], in 
which the MCDS1H and MC1H approximations have 
been implemented in the same way than for the MC- 
srDFT method [70, 71, 93], using the direct restricted- 
step second-order MCSCF algorithm of Jensen and 
coworkers [94-99]. For E x [n] and £7 c [n], we use two 
GGA exchange-correlation density functionals, Perdew- 
Burkc-Ernzcrhof (PBE) [100] and Becke/Lee- Yang-Parr 
(BLYP) [101, 102], without spin-density dependence. For 



implementing the density-scaled correlation functionals 
in the MCSCF algorithm, we need the scaling relations 
for the energy density, and its first- and second-order 
derivatives that we give in Appendix A. The computa- 
tional cost of the method is essentially the same as for 
a standard MCSCF calculation, with a small extra cost 
due to the DFT contribution. 

A good value for the empirical parameter A in Eqs. (7) 
and (8) is determined on the 03ADD6 benchmark 
set [103, 104] for the 1,3-dipolar cycloaddition reac- 
tions of ozone (O3) with ethylene (C2H4) or acetylene 
(C2H2) [90, 105, 106]. For these two reactions, there 
are three stationary points along the reaction coordi- 
nate: the van der Waals complex, the transition state, 
and the cycloadduct (primary ozonidc), all in a closed- 
shell spin-singlet state. The 03ADD6 set consists of the 
six energies of these stationary points of the two reac- 
tions, calculated relative to the energy of the separated 
reactants, and without zero-point vibrational energy cor- 
rection. Accurate calculations of these energies are dif- 
ficult and require to handle the subtle balance between 
static and dynamic correlation effects along the reaction 
coordinate. The ozone reactant, the van der Waals com- 
plex, and the transition state have a strong multicon- 
figurational character corresponding to the HOMO — > 
LUMO double excitation in ozone. In the cycloadduct, 
and to a less extent in the transition state, the stabi- 
lization of the ozone HOMO and the destabilization of 
the ozone LUMO greatly reduce this multiconfigurational 
character. In addition, there are small near-degeneracy 
correlation effects due to the ir and ir* orbitals of the 
reactive tt bond of ethylene and acetylene. For each sep- 
arate reactant, a CAS(2,2) wave function is chosen, the 
active space corresponding to the HOMO and LUMO or- 
bitals for ozone, and to the HOMO (tt) and LUMO (ir*) 
orbitals of the reactive tt bond for ethylene and acety- 
lene. For the van der Waals complex, the transition state, 
and the cycloadduct, a CAS(4,4) wave function is con- 
sistently chosen, the active space corresponding to the 
orbitals that connect to the ones chosen for the reactants 
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TABLE I: Energies of the van der Waals (vdW) complex, the transition state (TS), and the cycloadduct (primary ozonide), 
relative to the separated reactants, and the corresponding MAEs (in kcal/mol) for the addition of ozone with acetylene or 
ethylene (03ADD6 set), calculated by several methods. For the DS1DH-PBE and 1DH-BLYP double-hybrid approximations, 
we use the value A = 0.65 which was previously optimized in Ref. 79. For the MCDS1H and MC1H multiconfigurational hybrid 
approximations, we use a value of A = 0.25 which roughly minimizes the MAE according to Fig. 1. For the range-separated 
MC-srPBE multiconfigurational hybrid approximation, we use the value of the range-separation parameter fi — 0.40 bohr -1 
which was previously determined in Ref. 71. For the multiconfigurational methods, a complete active space wave function with 
2 electrons in 2 orbitals [CAS (4,4)] is chosen for the van der Waals complex, the transition state, and the cycloadduct, and 
a CAS(2,2) wave function for each separate reactant. All calculations were carried out with the aug-cc-pVTZ basis set. All 
calculations were done for M05/6-311+G(2df,2p) geometries, except for the MRMP2 results which are for CCSD(T)/cc-pVTZ 
geometries. 



3 + C2H2 -)• 3 + C2H4 -> 



Method 


vdW 


TS 


cycloadduct 


vdW 


TS 


cycloadduct 


MAE 


HF 


0.68 


23.08 


-87.12 


1.90 


17.91 


-82.58 


14.18 


MCSCF 


0.69 


27.54 


-77.25 


1.32 


22.13 


-68.06 


11.46 


MP2 


-3.18 


1.13 


-54.81 


-4.01 


-5.74 


-51.18 


5.67 


MRMP2 a 


-2.16 


8.77 


-48.19 


-2.09 


3.43 


-43.32 


5.16 


PBE 


-1.71 


-1.66 


-62.44 


-2.50 


-4.77 


-51.45 


4.23 




-0.57 


2.21 


-53.98 


-1.29 


-1.55 


-43.19 


6.03 


Single-hybrid approximations 
















PBEO 


-1.26 


1.65 


-74.00 


-1.55 


-1.74 


-64.74 


5.00 


B1LYP 


-0.60 


4.78 


-66.56 


-0.83 


0.71 


-57.47 


1.85 


B3LYP 


-0.67 


3.81 


-65.10 


-1.01 


-0.12 


-55.64 


2.06 


Double-hybrid approximations 
















DS1DH-PBE 


-2.08 


3.47 


-61.21 


-2.54 


-1.81 


-54.90 


2.51 


1DH-BLYP 


-1.92 


4.42 


-58.37 


-2.37 


-1.14 


-52.12 


3.12 


B2-PLYP 6 


-1.47 


5.00 


-60.18 


-1.81 


-0.13 


-53.05 


2.42 


Multiconfigurational hybrid approximations 
















MC-srPBE 


-0.93 


4.12 


-72.73 


-0.87 


0.71 


-65.53 


4.27 


MCDS1H-PBE 


-1.06 


3.88 


-70.41 


-1.22 


0.30 


-60.71 


3.11 


MC1H-PBE 


-1.08 


3.66 


-70.97 


-1.25 


0.13 


-61.26 


3.35 


MCDS1H-BLYP 


0.28 


7.94 


-62.47 


0.26 


3.78 


-52.86 


1.77 


MC1H-BLYP 


-0.36 


6.74 


-63.76 


-0.47 


2.57 


-54.21 


1.30 


Best estimate" 


-1.90 


7.74 


-63.80 


-1.94 


3.37 


-57.15 





"From Ref. 90. 



^Performed with the Gaussian09 program [91]. 



in the dissociation limit. We use the aug-cc-pVTZ basis 
set [107, 108] and the fixed geometries of Ref. 90 opti- 
mized using the hybrid meta-GGA exchange-correlation 
functional M05 [109] with the 6-311+G(2df,2p) basis 
set [110, 111]. The reference values for the energies are 
from Ref. 90 and were obtained from extensive coupled- 
cluster calculations extrapolated to the complete basis 
set limit [90, 105]. We calculate the mean absolute error 
(MAE) over the six values as a function of the parameter 
A. Wc compare the MCDS1H and MC1H approxima- 
tions with (a) some non-hybrid methods: HF, MCSCF, 
MP2 [112], multireference MP 2 (MRMP2) [90, 113], 
PBE [100], and BLYP [101, 102]; (b) some single- 
hybrid approximations: PBEO [114, 115], B1LYP [116], 
and B3LYP [117, 118]; (c) some double-hybrid approx- 
imations: DS1DH-PBE [79], 1DH-BLYP [79], and B2- 
PLYP [89], all applied in a spin- restricted formalism. 
We also compare with the range-separated MC-srPBE 



multiconfigurational hybrid approximation [70-72] using 
the short-range PBE exchange-correlation functional of 
Ref. 119 and the value of the range-separation parame- 
ter fi = 0.40 bohr -1 which was previously determined in 
Ref. 71. 

We also test the MCDS1H and MC1H approximations 
by computing the potential energy curves the five di- 
atomic molecules H2, Li2, C2, N2, and F2, using in each 
case a full-valence CAS wave function and the cc-pVTZ 
basis set [107]. 



IV. RESULTS 

A. 03ADD6 database 

Figure 1 shows the MAEs for the 03ADD6 set as func- 
tions of the parameter A for the MCDS1H and MC1H 
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FIG. 2: Potential energy curves of H2 calculated by HF, MCSCF, and several methods based on the PBE (left) or BLYP 
(right) exchange-correlation density functionals. For the MCDS1H and MC1H multiconfigurational hybrid approximations, we 
use a value of A = 0.25. For the range-separated MC-srPBE multiconfigurational hybrid approximation, we use a value of the 
range-separation parameter of fi = 0.40 bohr -1 . For all multiconfigurational methods, we use a full- valence CAS wave function. 
The basis set is cc-pVTZ basis. The accurate curve is from Ref. 28. 



approximations with the BLYP and PBE exchange- 
correlation density functionals. For A = 0, the MCDS1H 
and MC1H approximation reduce to a standard KS calcu- 
lation with the corresponding approximate density func- 
tional. For A = 1, they reduce to a standard MCSCF 
calculation. Contrary to what was found for the calcu- 
lation of atomization energies using double-hybrid ap- 
proximations [79], here neglecting density scaling in the 
correlation functional makes little difference. Toward the 
A = 1 end of the curves, the MCDS1H and MC1H ap- 
proximations inherit the inaccuracy of MCSCF which 
neglects dynamic correlation. The MAE curves of the 
MCDS1H-BLYP and MC1H-BLYP approximations dis- 
play a marked minimum at an intermediate value of A, 
thus improving upon both the standard BLYP and MC- 
SCF calculations. The minimum is reached at A = 0.25 
for MCDS1H-BLYP and at A = 0.30 for MC1H-BLYP, 
with MAEs below 2 kcal/mol. For the MCDS1H-PBE 
and MC1H-PBE approximations, the MAE curves have 
a plateau around A = 0.25 with a MAE of about 3 
kcal/mol, which is again smaller than both the standard 
PBE and MCSCF calculations. In view of these results, 
we choose the value A = 0.25 in all MCDS1H and MC1H 
approximations. It is a conservative choice since it gives 
the same fraction of exact exchange as the one usually ad- 
vocated in the usual single-hybrid approximations [120]. 

Table I reports the energies of the van der Waals com- 
plex, the transition state, and the cycloadduct of the two 
reactions of the 03ADD6 set, relative to the separated 
reactants, calculated with the MCDS1H and MC1H ap- 
proximations at A = 0.25. For comparison, we also report 
results for various non-hybrid and other hybrid methods. 

During the early stages of the two reactions, a weakly 
bound van der Waals complex is formed which lies in 
a shallow minimum (-1.90 kcal/mol for C2H2 and -1.94 
kcal/mol for C 2 H 4 ) below the reactants. The MCDS1H 
and MC1H approximations give significantly underesti- 
mated well depths, which are still in improvement over 



standard MCSCF but not over standard KS calculations 
with the corresponding functionals. The range-separated 
MC-srPBE method does also not perform better than KS 
PBE for these van der Waals systems. A better descrip- 
tion of the long-range dispersion correlation would indeed 
require inclusion of perturbation corrections on top of 
the active space [73]. As expected, the double- hybrid ap- 
proximations, which include second-order perturbation 
corrections, tend to perform better for these van der 
Waals complexes. The best performance is achieved with 
MRMP2 which is able to correctly describe both multi- 
configurational effects and dispersion correlations. How- 
ever, one should keep in mind than the methods MP2, 
MRMP2, and the double hybrids are most likely less con- 
verged with respect to the basis size than the other meth- 
ods. 

The activation barriers of the transition states are un- 
derestimated (or not present at all) in KS PBE and BLYP 
calculations, and to a less extent with the single-hybrid 
and double-hybrid approximations, while they are largely 
overestimated in MCSCF. Note that here, contrary to 
the common case, MCSCF gives higher activation barri- 
ers than HF because the ozone reactant has more static 
correlation than the transition state. The MCDS1H and 
MC1H approximations give an improvement of about 5 
kcal/mol over the KS calculations with the correspond- 
ing functionals. The range-separated MC-srPBE method 
gives activation barriers which are slightly better than 
the ones given by MCDS1H-PBE and MC1H-PBE, but 
largely worse than the ones given by MCDS1H-BLYP 
and MC1H-BLYP. The values obtained with MCDS1H- 
BLYP and MC1H-BLYP, as well as with MRMP2, are 
all within 1 kcal/mol of the best estimates. 

The reaction energies of the formation of the cy- 
cloadducts are overestimated in MCSCF and underes- 
timated in MRMP2 (by about 15 kcal/mol). Zhao et 
al. observed that even using a large (14,14) active space 
does not improve the MRMP2 reaction energy [90]. All 
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FIG. 3: Potential energy curves of L12, C2, N2, and F2 calculated with MCSCF and several methods based on the PBE 
(left) or BLYP (right) exchange-correlation density functionals. For the MCDSfH and MC1H multiconfigurational hybrid 
approximations, we use a value of A = 0.25. For the range-separated MC-srPBE multiconfigurational hybrid approximation, we 
use a value of the range-separation parameter of /i = 0.40 bohr -1 . For all multiconfigurational methods, we use a full- valence 
CAS wave function. The basis set is cc-pVTZ basis. The accurate curves are from Ref. 28. 
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the hybrid methods give more reasonable reaction ener- 
gies. In particular, MC1H-BLYP gives reaction energies 
within less than 3 kcal/mol of the best estimates. 

If we accept to look more closely to the MAE values in 
spite of the limited statistics, we see that MC1H-BLYP 
gives overall the smallest MAE with 1.30 kcal/mol. The 
multiconfigurational hybrid MCDS1H-PBE involving the 
PBE exchange-correlation functional gives a larger MAE 
of 3.11 kcal/mol, but turns out to perform better on aver- 
age than the range-separated MC-srPBE method which 
gives a MAE of 4.27 kcal/mol. Neglecting density scal- 
ing in the correlation functional of multiconfigurational 
hybrids seems slightly favorable for BLYP and slightly 
unfavorable for PBE. This is in line with what was found 
for double-hybrid approximations [79] , although the MP2 
correlation part made it much more sensitive to the ne- 
glect of the density scaling. Even if the effect of neglect- 
ing density scaling is systematic in giving more nega- 
tive complement correlation energies [n] for all density 
functionals, its effect on the MAE depends more on for- 
tuitous compensation of errors for the approximate func- 
tional used. 



B. Dissociation of diatomic molecules 

We now turn to the calculation of potential energy 
curves of diatomic molecules. This is a harder problem 
since static correlation effects are dominant at dissoci- 
ation. For the multiconfigurational hybrids, we report 
here only the curves of MCDS1H for PBE and of MC1H 
for BLYP according to the results of Section IV A, but the 
differences between the curves of MCDS1H and MC1H 
for both BLYP and PBE arc in fact very small for these 
diatomic molecules. 

Figure 2 shows the potential energy curve of H2 calcu- 
lated by hybrid approximations using the PBE and BLYP 
density functionals. Around the equilibrium internuclear 
distance, all DFT-based method, including MCDS1H- 
PBE and MC1H-BLYP, are accurate, which means that 
they properly describe dynamic correlation. At large dis- 
tances, the single-hybrid approximations (PBEO, B1LYP, 
and B3LYP), which include a fraction of HF exchange 
energy, give less accurate potential energy curves than 
non-hybrid KS calculations (PBE and BLYP). By in- 
clusion of a fraction of exact static correlation energy, 
the multiconfigurational hybrids (MCDS1H-PBE, MC- 
srPBE, and MC1H-BLYP) correct this behavior and give 
potential energy curves that correctly saturate beyond a 
distance of about 5 bohr, as the MCSCF curve does. This 
point is explained by a detailed analysis of the asymp- 
totic expansion of the potential energy curves in a mini- 
mal basis in Appendix B. However, the MCDS1H-PBE, 
MC-srPBE, and MC1H-BLYP methods still display sig- 
nificant errors on the energy of the separated atoms due 
to the density-functional approximations. Indeed, as in 
restricted KS calculations, the density functionals used in 
the multiconfigurational hybrids depend only of the to- 



tal density and do not give accurate energies in the limit 
of separated atoms of open-shell character. In an unre- 
stricted KS calculation, the energy at dissociation can 
be improved by breaking the spin symmetry and there- 
fore introducing a fictitious spin density which helps to 
describe the separated atoms. In our present implemen- 
tation of the multiconfigurational hybrids, the spin sym- 
metry is imposed on the MCSCF wave function so that 
there is no fictitious spin density to be used in the density 
functionals. 

Figure 3 shows the potential energy curves of Li2, C2, 
N2, and F2. The results are similar than for H2. Around 
the equilibrium distance, the MCDS1H-PBE and MC1H- 
BLYP curves are similar to the standard hybrid or non- 
hybrid KS calculations. At large distances, the single- 
hybrid approximations give a spurious increase of the 
energy, whereas the MCDS1H-PBE and MC1H-BLYP 
curves correctly saturate. The MC1H-BLYP approxi- 
mation gives good total energies, but again a signifi- 
cant error remains at dissociation, especially for N2 . The 
DS1DH-PBE and MC-srPBE approximation gives curves 
of very similar shape. 



V. CONCLUSIONS 

We have presented a multiconfigurational hybrid 
density-functional theory which rigorously combines MC- 
SCF and DFT based on a linear decomposition of the 
electron-electron interaction. It is straightforward exten- 
sion of the usual hybrid approximations by essentially 
adding a fraction A of exact static correlation in addi- 
tion to the fraction A of exact exchange. Any existing 
approximate exchange-correlation density functional can 
be used in this scheme by using a simple scaling relation 
with A. Test calculations on the cycloaddition reactions 
of ozone with ethylene or acetylene and the dissociation 
of diatomic molecules with the PBE and BLYP density 
functionals show that a good value of A is 0.25, as in the 
usual hybrid approximations. 

Interestingly, the results seem to indicate that the 
present approach based on a simple linear decomposition 
of the electron-electron interaction is at least as good as 
the range-separated multiconfigurational hybrid method 
of Ref. 70-72 for including static correlation in DFT, at 
least with the approximate density functionals used here. 
Of course, with better short-range density-functional ap- 
proximations (in particular, we do not have a short-range 
version of the LYP correlation functional for compari- 
son), the conclusion could be different. Also, one should 
note that hybrid approaches combining perturbation the- 
ory with DFT based on a linear decomposition of the 
interaction [79] do not have the advantages of the range- 
separated hybrid approaches for fast basis-size conver- 
gence and explicit inclusion of long-range van der Waals 
interactions. For MCSCF, however, basis set convergence 
is not so much an issue. 

The present results suggest that the proposed multi- 
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configurational hybrid approximations can improve over 
usual DFT approximations for situations with strong 
static correlation effects. It remains however to assess the 
performance of this multiconfigurational hybrid method 
on a larger variety of systems. Future work includes 
adding the dependence on the spin density in the func- 
tionals to be able to properly handle open-shell systems, 
and possibly other additional variables such as the on- 
top pair density as an alternative to the spin density for 
improving the accuracy of closed-shell systems [121]. 
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Appendix A: Scaling relations for the derivatives of 
the density-scaled correlation functional 

We give the scaling relations for the density-scaled cor- 
relation functional ££[n] = \ 2 E c [n 1 /x] and its deriva- 
tives in the case of generalized-gradient approximations 
(GGA). Starting from a standard GGA density func- 
tional written as 

£ c ,ggaM = /e c (n(r),|Vn(r)|)dr, (Al) 

where |Vn(r)| is the norm of the density gradient, the 
corresponding scaled functional is 

£c,ggaN = J e^(n(r),\Vn(r)\)dr, (A2) 

where the energy density is obtained by scaling relation 
(see Rcf. 79) 

^(„(r) J |Vn(r)|) = A«e f!#,J^WlV (A3) 



A 3 A 4 J 
The first-order derivatives of the energy density are 

M Wr) , |Vn(r)l) .^(^),M, (A4) 



and 

de X c 
0|Vn| 



(n(r),|Vn(r)|) = A- 



de c ( n(r) |Vn(r)| 



9|Vn| V A 3 ' A 4 



The second-order derivatives are 



32 „ A 



dn 2 



(n(r),|Vn(r)|) = -- 



1 d 2 e c fn(r) |Vn(r) 



A dn 2 V A 3 ' A 4 



(A5) 



,(A6) 



8 2 e x c 
d\Vn\ 
1 



(n(r),|Vn(r)|) 



A 3 a|Vn| 2 



n(r) |Vra(r)| 



A 3 



A 4 



(A7) 



and 



2 „ A 



d 2 e 



dnd | Vn 
1 d 2 



(n(r),|Vn(r)|) = 

n(r) |Vn(r) 



A 2 dnd\Vn\ 



A 3 



A 4 



(A8) 



Appendix B: Asymptotic expansion of the potential 
energy curve of H2 

We consider the H2 molecule in a Slater minimal basis, 
with a basis function a localized on the left atom and a 
basis function b localized on the right atom, both basis 
functions being identical with exponent £ = 1. In the 
large internuclear distance R limit, the two molecular 
orbitals are 1 = (a + b)/y/2 and 2 = (a - b)/y/2. The 
total restricted Hartree-Fock (RHF) energy writes 



E 



RHF 



= 2h n + Jn 



1 



(Bl) 



where h\\ = t\\ + V\\ is the sum of the kinetic inte- 
gral in = (l|i|l) and the nuclei-electron integral «u = 
(l|w„ e |l), and Jn = (H|H) is the Coulomb two-electron 
integral, and 1/R is the nuclear repulsion energy. By ex- 
panding the molecular orbital 1 into the localized func- 
tions a and b, and using the symmetry between a and b, 
it is easy to find the large R behavior of all these terms: 



(a\t\a) + (a\i\b) 



+ 0(e- R ), 



(B2) 



and 



vu = (a\v ne \a) + (a\v ne \b) = -1 - - + 0{e- H ), (B3) 

H 



and 



(aa\aa) (aa\bb) 
~^~ + ^^~ 

- + — + 0(e- R ) 
16 2R v ' 



2(aa\ab) + (ab\ab) 



(B4) 



where 0(e~ R ) stands for exponentially decaying terms in 
R. For the values of the integrals, see Ref. 122. Adding 
all the pieces together, it leads to the following asymp- 
totic expansion of the total RHF energy 



£RHF = _ I! _ J_ + Q , e - R) 

16 2R V ' 



(B5) 



At dissociation, the RHF wave function contains 50% 
of the incorrect ionic contribution H + ...H~, which is re- 
sponsible for too high an energy and for the spurious 
electrostatic attraction term —1/2R. 
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The full configuration interaction (FCI) correlation en- 
ergy in this basis is found by diagonalizing the 2x2 
Hamiltonian matrix, leading to 



1 



2 ' 



(B6) 

where E2 = 2h 22 + J22 + 1/-R is the energy of the double- 
excited determinant, and Ki 2 = (12 1 12) is the exchange 
two-electron integral. The asymptotic behavior of E 2 is 
exactly the same as the one of £; RHF , so that E 2 — £? RHF 
vanishes exponentially when R — > 00 and the asymptotic 
behavior of E FCl is determined by Ki 2 only: E FCl = 
—K\ 2 + 0(e~ R ). The asymptotic behavior of K\ 2 is 



(aa\aa) (aa\bb) 5 1 _ R 
Al2 = — — = -Q-2R + ° {e ] 



giving for the correlation energy 



Ej cl = -- + — +0(e- R ). 

16 2R v ' 



(B7) 



(B8) 



Adding the asymptotic expansions of Eqs. (B5) and (B8) 
gives the asymptotic expansion of the total FCI energy 
in this basis 



E FCI = -1 + 0(e- H ) 



-R\ 



(B9) 



which implies that the FCI potential energy curve satu- 
rates quickly at large internuclear distance. 

In restricted Kohn-Sham (RKS) density-functional 
theory, the total energy writes 



E 



RKS 



2hn + 2J\\ + E x + E c 



1 



(BIO) 



where E x and E c are the exchange and correlation ener- 
gies. At large R, it behaves as 



E 



RKS 



+ E x +E c + 0{e- R ). 



(Bll) 



With local or semilocal density-functional approxima- 
tions, E x and E c go exponentially to constants when 
R —> 00, so that the asymptotic expansion of _E RKS does 
not contain a spurious term in 1/R. 

Single-hybrid approximations introduces a fraction A 
of RHF exchange which have the following asymptotic 
expansion 



E 



RHF 



-Ju 



5 

16 



2R 



+ 0(e- R ), 



(B12) 



and therefore introduce a wrong —\/2R term in the total 
energy, 



^brid = _ 3 _ 5A + (i _ x)Ex + Ec _A_ + (e- R ). 

(B13) 

Single-hybrid approximations thus deteriorate the large 
R behavior of local or semilocal density-functional ap- 
proximations (see Fig. 2). The multiconfigurational hy- 
brid approximations introduced in this work correct this 
behavior by adding a fraction of the FCI correlation en- 
ergy which, in the limit of large R, is just XE FCl , the 
linearity in A being a signature of static correlation. For 
example, the MC1H approximation has the following 
asymptotic expansion 



E 



31H = - | - y + (1 - X)E X + (1 - \ 2 )E C + 0(e~ R ), 



(B14) 

with no longer any spurious 1/R term, and thus improves 
the large R behavior (see Fig. 2). It is a typical exam- 
ple where exact exchange and static correlation must be 
considered together. 

For range-separated density-functional theory, the sit- 
uation is similar. Range-separated single-hybrid approx- 
imations [80, 123-126] include some long-range RHF ex- 
change and their asymptotic expansions display a wrong 
—1/2R term, just as RHF. Their behavior for large R 
is in fact worse than that of usual single-hybrid approx- 
imations since the —1/2R term is not weighted by A. 
However, the range-separated CI-srDFT [66, 67] or MC- 
srDFT [70-72] methods add some exact long-range cor- 
relation energy which removes this wrong —1/2R term. 

Note that other forms of single-hybrid approximations 
which do not use RHF exchange at long range [127, 128] 
allow one to avoid a wrong —1/2R term in the large 
R limit. Symmetry breaking is another way to avoid 
a wrong asymptotic — 1/2_R term since the unrestricted 
Hartree-Fock exchange energy does not contain such a 
term. 

The fact that local or semilocal approximations for E x 
and E c do not introduce 1/R terms is in agreement with 
the usual conviction that approximate CCA exchange 
functionals not only represent exchange but also static 
correlation, while approximate GGA correlation func- 
tionals represent dynamic correlation only [4] . 



[1] P. Hohenberg and W. Kohn, Phys. Rev. 136, B 864 [3] W. Koch and M. C. Holthausen, A Chemist's Guide 

(1964) . To Density Functional Theory (Wiley- VCH, New York, 
[2] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 2001). 

(1965) . [4] O. V. Gritsenko, P. R. T. Schipper, and E. J. Baerends, 



10 



J. Chcm. Phys. 107, 5007 (1997). 
[5] A. D. Becke, J. Chem. Phys. 119, 2972 (2003). 
[6] A. J. Cohen, P. Mori-Sanchez, and W. Yang, Science 

321, 792 (2008). 
[7] N. E. Schultz, Y. Zhao, and D. G. Truhlar, J. Phys. 

Chem. A 109, 11127 (2005). 
[8] D. Cremer, Mol. Phys. 99, 1899 (2001). 
[9] J. Slater, J. Mann, T. Wilson, and J. Wood, Phys. Rev. 
184, 672 (1969). 
[10] M. Levy, Phys. Rev. A 26, 1200 (1982). 
[11] E. H. Lieb, Int. J. Quantum. Chem. 24, 24 (1983). 
[12] B. I. Dunlap and W. N. Mei, J. Chem. Phys. 78, 4997 
(1983). 

[13] S. G. Wang and W. H. E. Schwarz, J. Chem. Phys. 105, 
4641 (1996). 

[14] P. R. T. Schipper, O. V. Gritsenko, and E. J. Baerends, 

Theor. Chem. Acc. 99, 329 (1998). 
[15] J. D. Goddard and G. Orlova, J. Chem. Phys. Ill, 7705 

(1999). 

[16] M. Filatov and S. Shaik, Chem. Phys. Lett. 304, 429 
(1999). 

[17] M. Filatov and S. Shaik, J. Phys. Chem. A 103, 8885 
(1999). 

[18] M. Filatov and S. Shaik, Chem. Phys. Lett. 332 , 409 

(1999) . 

[19] M. Filatov, S. Shaik, M. Woeller, S. Grimme, and 

S. Peyerimhoff, Chem. Phys. Lett. 316, 135 (2000). 
[20] M. Filatov and S. Shaik, J. Phys. Chem. A 104, 6628 

(2000) . 

[21] C. A. Ullrich and W. Kohn, Phys. Rev. Lett. 87, 093001 

(2001) . 

[22] K. Y. R. Takeda, S. Yamanaka, Int. J. Quantum. Chem. 

93, 317 (2003). 
[23] S. Grimme and M. Waletzke, J. Chem. Phys. Ill, 5645 

(1999). 

[24] E. V. Beck, E. A. Stahlberg, L. W. Burggraf, and J.-P. 

Blaudeau, Chem. Phys. 349, 158 (2008). 
[25] Q. Wu, C.-L. Cheng, and T. Van Voorhis, J. Chem. 

Phys. 127, 164119 (2007). 
[26] Q. Wu, B. Kaduk, and T. Van Voorhis, J. Chem. Phys. 

130, 034109 (2009). 
[27] G. C. Lie and E. Clementi, J. Chem. Phys. 60, 1275 

(1974). 

[28] G. C. Lie and E. Clementi, J. Chem. Phys. 60, 1288 

(1974) . 

[29] R. Colle and O. Salvetti, Theor. Chim. Acta 37, 329 

(1975) . 

[30] R. Colle and O. Salvetti, Theor. Chim. Acta 53, 55 
(1979). 

[31] R. Colle and O. Salvetti, J. Chem. Phys. 79, 1404 
(1983). 

[32] R. Colle and O. Salvetti, J. Chem. Phys. 93 (1990). 
[33] A. Savin, Int. J. Quantum. Chem. 22, 59 (1988). 
[34] A. Savin, J. chim. phys. 86, 757 (1989). 
[35] A. Savin, in Density functional methods in chemistry, 

edited by J. Labanowski and J. Andzelm (Springer- 

Verlag, New York, 1991), p. 213. 
[36] B. Miehlich, H. Stoll, and A. Savin, Mol. Phys. 91, 527 

(1997). 

[37] C. Gutle and A. Savin, Phys. Rev. A 75, 032519 (2007). 
[38] F. Moscardo and E. San-Fabian, Int. J. Quantum. 

Chem. 40, 23 (1991). 
[39] E. Kraka, Chem. Phys. 161, 141 (1992). 
[40] E. Kraka, D. Cremer, and S. Nordholm, in Molecules 



in Natural Science and Biomedicme, edited by Z. Mak- 
sic and M. Eckert-Maksic (Ellis Horwood, Chichester, 
1992), p. 351. 

[41] W. Wu and S. Shaik, Chem. Phys. Lett. 301, 37 (1999). 

[42] H. Stoll, Chem. Phys. Lett. 376, 141 (2003). 

[43] F. Ying, P. Su, Z. Chen, S. Shaik, and W. Wu, J. Chem. 

Theory Comput. 8, 1608 (2012). 
[44] N. O. J. Malcolm and J. J. W. McDouall, J. Phys. 

Chem. 98, 12579 (1994). 
[45] N. O. J. Malcolm and J. J. W. McDouall, J. Phys. 

Chem. 100, 10131 (1996). 
[46] N. O. J. Malcolm and J. J. W. McDouall, J. Phys. 

Chem. A 101, 8119 (1997). 
[47] N. O. J. Malcolm and J. J. W. McDouall, Chem. Phys. 

Lett. 282, 121 (1998). 
[48] J. J. W. McDouall, Mol. Phys. 101, 361 (2003). 
[49] P. Borowski, K. D. Jordan, J. Nichols, and P. Nachtigall, 

Theor. Chem. Acc. 99, 135 (1998). 
[50] J. Grafenstein and D. Cremer, Chem. Phys. Lett. 316, 

569 (2000). 

[51] J. Grafenstein and D. Cremer, Phys. Chem. Chem. 

Phys. 2, 2091 (2000). 
[52] J. Grafenstein and D. Cremer, Mol. Phys. 103, 279 

(2005) . 

[53] R. Takeda, S. Yamanaka, and K. Yamaguchi, Chem. 

Phys. Lett. 366, 321 (2002). 
[54] K. Nakata, T. Ukai, S. Yamanaka, T. Takada, and 

K. Yamaguchi, Int. J. Quantum. Chem. 106, 3325 

(2006) . 

[55] S. Yamanaka, K. Nakata, T. Ukai, T. Takada, and 
K. Yamaguchi, Int. J. Quantum. Chem. 106, 3312 
(2006). 

[56] T. Ukai, K. Nakata, S. Yamanaka, T. Takada, and 

K. Yamaguchi, Mol. Phys. 105, 2667 (2007). 
[57] S. Gusarov, P.-A. Malmqvist, R. Lindh, and B. O. Roos, 

Theor. Chim. Acc. 112, 84 (2004). 
[58] A. J. Perez-Jimenez and J. M. Perez- Jorda, Phys. Rev. 

A 75, 012503 (2007). 
[59] A. J. Perez- Jimenez, J. M. Perez- Jorda, I. D. P. R. Mor- 

eira, and F. Illas, J. Comput. Chem. 28, 2559 (2007). 
[60] A. J. Perez- Jimenez, J. M. Perez- Jorda, and J. C. 

Sancho-Garda, J. Chem. Phys. 127, 104102 (2007). 
[61] M. Weimer, F. D. Sala, and A. Gorling, J. Chem. Phys. 

128, 144109 (2008). 
[62] Y. Kurzweil, K. V. Lawler, and M. Head-Gordon, Mol. 

Phys. 107, 2103 (2009). 
[63] A. Savin and H.-J. Flad, Int. J. Quantum. Chcm. 56, 

327 (1995). 

[64] A. Savin, in Recent Advances in Density Functional 

Theory, edited by D. P. Chong (World Scientific, 1996). 
[65] A. Savin, in Recent Developments of Modern Density 

Functional Theory, edited by J. M. Seminario (Elsevier, 

Amsterdam, 1996), pp. 327-357. 
[66] T. Leininger, H. Stoll, H.-J. Werner, and A. Savin, 

Chem. Phys. Lett. 275, 151 (1997). 
[67] R. Pollet, A. Savin, T. Leininger, and H. Stoll, J. Chem. 

Phys. 116, 1250 (2002). 
[68] A. Savin, F. Colonna, and R. Pollet, Int. J. Quantum. 

Chem. 93, 166 (2003). 
[69] J. Toulouse, F. Colonna, and A. Savin, Phys. Rev. A 

70, 062505 (2004). 
[70] J. K. Pedersen and H. J. A. Jensen, a second order 

MCSCF-DFT hybrid algorithm (unpublished). 
[71] E. Fromager, J. Toulouse, and H. J. A. Jensen, J. Chem. 



11 



Phys. 126, 074111 (2007). 
[72] E. Fromager, F. Real, P. Wahlin, U. Wahlgren, and 

H. J. A. Jensen, J. Chem. Phys. 131, 054107 (2009). 
[73] E. Fromager, R. Cimiraglia, and H. J. A. Jensen, Phys. 

Rev. A 81, 024502 (2010). 
[74] T. Tsuchimochi, G. E. Scuseria, and A. Savin, J. Chem. 

Phys. 132, 024111 (2010). 
[75] T. Tsuchimochi and G. E. Scuseria, J. Chem. Phys. 134, 

064101 (2011). 
[76] K. Pernal, Phys. Rev. A 81, 052511 (2010). 
[77] D. R. Rohr, J. Toulouse, and K. Pernal, Phys. Rev. A 

82, 052502 (2010). 
[78] D. R. Rohr and K. Pernal, J. Chem. Phys. 135, 074104 

(2011). 

[79] K. Sharkas, J. Toulouse, and A. Savin, J. Chem. Phys. 

134, 064113 (2011). 
[80] J. G. Angyan, I. C. Gerber, A. Savin, and J. Toulouse, 

Phys. Rev. A 72, 012510 (2005). 
[81] M. Levy and J. P. Perdew, Phys. Rev. A 32, 2010 

(1985) . 

[82] M. Levy, W. Yang, and R. G. Parr, J. Chem. Phys. 83, 

2334 (1985). 
[83] M. Levy, Phys. Rev. A 43, 4637 (1991). 
[84] M. Levy and J. P. Perdew, Phys. Rev. B 48, 11638 

(1993). 

[85] A. D. Becke, J. Chem. Phys. 104, 1040 (1996). 

[86] M. Ernzerhof, J. P. Perdew, and K. Burke, in Density 

Functional Theory, edited by R. Nalewajski (Springer- 

Verlag, Berlin, 1996). 
[87] E. Fromager and H. J. A. Jensen, Phys. Rev. A 78, 

022504 (2008). 
[88] J. G. Angyan, Phys. Rev. A 78, 022510 (2008). 
[89] S. Grimme, J. Chem. Phys. 124, 034108 (2006). 
[90] Y. Zhao, O. Tishchenko, J. R. Gour, W. Li, J. J. Lutz, 

P. Piecuch, and D. G. Truhlar, J. Phys. Chem. A 113, 

5786 (2009). 

[91] M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. 
Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, 
V. Barone, B. Mennucci, G. A. Petersson, et al., Gaus- 
sian 09 Revision A.l, Gaussian Inc. Wallingford CT 
2009. 

[92] C. Angeli, K. L. Bak, V. Bakken, O. Christiansen, 
R. Cimiraglia, S. Coriani, P. Dahle, E. K. Dalskov, 
T. Enevoldsen, B. Fernandez, et al., Dalton, a molecular 
electronic structure program, release dalton2011 (2011), 
see http://daltonprogram.org. 

[93] J. K. Pedersen, PhD thesis, University of Southern Den- 
mark, Odense (2004). 

[94] H. J. A. Jensen and P. J0rgensen, J. Chem. Phys. 80, 
1204 (1984). 

[95] H. J. A. Jensen and H. Agren, Chem. Phys. Lett. 110, 
140 (1984). 

[96] H. J. A. Jensen and H. Agren, Chem. Phys. 104, 229 

(1986) . 

[97] H. J. A. Jensen, P. Jorgensen, and H. Agren, J. Chem. 

Phys. 87, 451 (1987). 
[98] H. J. A. Jensen, P. J0rgensen, H. Agren, and J. Olsen, 

J. Chem. Phys. 88, 3834 (1988). 
[99] H. J. A. Jensen, in Relativistic and Electron Correlation 



Effects in Molecules and Solids, edited by G. L. Malli 

(Plenum, New York, 1994), p. 179. 
[100] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 
101] A. D. Becke, Phys. Rev. A 38, 3098 (1988). 
102] C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 37, 785 

(1988). 

103] L. Goerigk and S. Grimme, J. Chem. Theory Comput 

6, 107 (2010). 

104] L. Goerigk and S. Grimme, J. Chem. Theory Comput 

7, 291 (2011). 

105] S. E. Wheeler, D. H. Ess, and K. N. Houk, J. Phys. 

Chem. A 112, 1798 (2007). 
106] T. Saito, S. Nishihara, Y. Kataoka, Y. Nakanishi, 

Y. Kitagawa, T. Kawakami, SYamanaka, M. Oku- 

mura, and K. Yamaguchi, J. Phys. Chem. A 114, 12116 

(2010). 

107] T. H. Dunning, J. Chem. Phys. 90, 1007 (1989). 

108] R. A. Kendall, T. H. Dunning, and R. J. Harrison, J. 

Chem. Phys. 96, 6796 (1992). 
109] Y. Zhao, N. E. Schultz, and D. G. Truhlar, J. Chem. 

Phys. 123, 161103 (2005). 
110] R. Krishnan, J. S. Binkley, R. Seeger, and J. A. Pople, 

J. Chem. Phys. 72, 650 (1980). 
Ill] A. D. McLean and G. S. Chandler, J. Chem. Phys. 72, 

5639 (1980). 

112] C. M0ller and M. S. Plesset, Phys. Rev. 46, 618 (1934). 
113] K. Hirao, Chem. Phys. Lett. 190, 374 (1992). 
114] C. Adamo and V. Barone, J. Chem. Phys. 110, 6158 
(1999). 

115] C. Adamo and V. Barone, Chem. Phys. Lett. 298, 113 
(1998). 

116] C. Adamo and V. Barone, Chem. Phys. Lett. 274, 242 
(1997). 

117] A. D. Becke, J. Chem. Phys. 98, 5648 (1993). 
118] V. Barone and C. Adamo, Chem. Phys. Lett. 224, 432 
(1994). 

119] E. Goll, H.-J. Werner, H. Stoll, T. Leininger, P. Gori- 
Giorgi, and A. Savin, Chem. Phys. 329, 276 (2006). 

120] J. P. Perdew, M. Ernzerhof, and K. Burke, J. Chem. 
Phys 105, 9982 (1996). 

121] J. P. Perdew, A. Savin, and K. Burke, Phys. Rev. A 51, 
4531 (1995). 

122] M. J. S. Dewar and J. Kelemen, J. Chem. Educ. 48, 494 
(1971). 

123] H. Iikura, T. Tsuneda, T. Yanai, and K. Hirao, J. Chem. 

Phys. 115, 3540 (2001). 
124] T. Yanai, D. P. Tew, and N. C. Handy, Chem. Phys. 

Lett. 393, 51 (2004). 
125] I. C. Gerber and J. G. Angyan, Chem. Phys. Lett. 415, 

100 (2005). 

126] O. A. Vydrov and G. E. Scuseria, J. Chem. Phys. 125, 
234109 (2006). 

127] J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. 

Phys. 118, 8207 (2003). 
128] T. M. Henderson, A. F. Izmaylov, G. E. Scuseria, and 

A. Savin, J. Chem. Phys. 127, 221103 (2007). 



